r 3 



Europaisches Patentamt 
European Patent Office 
Office europeen des brevets 



(11) 



EP 1 203 560 A2 



(12) 



EUROPEAN PATENT APPLICATION 



(43) 


Date of nuhli cation* 


\p\) IntCI.'. MOID D/UO 




08 05 2002 Bulletin 2002/19 


(21) 


Application number: 01309190.5 




(22) 


Date of filing: 30.10.2001 




fQA\ 

(84) 


Designated Contracting States: 


(72) Inventor: Govan, Assaf 




AT BE CH CY DE DK ES Fl FR GB GR IE IT LI LU 


Haifa 34400 (IL) 




■ ■ KM III TO 

MC NL FT SE TR 






Designated Extension States: 


(74) Representative: Mercer, Christopher Paul et al 




AL LT LV MK RO SI 


Carpmaels & Ransford 






43, Bloomsbury Square 


(30) 


Priority: 31.10.2000 US 702973 


London WC1A2RA (GB) 


(71) 


Applicant: Biosense, Inc. 






New Brunswick, New Jersey 08933-7003 (US) 





(54) Method and apparatus for tracking an object, in particular a catheter, using energy fields 



(57) A method for tracking an object, including pro- 
ducing an unperturbed energy field at a plurality of pre- 
determined frequencies in the vicinity of the object and 
determining a characteristic of a perturbing energy field 
induced responsive to the unperturbed field, due to in- 
troduction of an article, responsive to the unperturbed 
field, into the vicinity of the object. The method further 
includes receiving a plurality of resultant signals respon- 



sive to the unperturbed and perturbing energy fields 
generated at a location of the object after introduction 
of the article, determining an optimal frequency for the 
unperturbed energy field from amongst the plurality of 
predetermined frequencies responsive to a parameter 
of the resultant signals, and determining spatial coordi- 
nates of the object responsive to the resultant signal at 
the optimal frequency. 
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Description 

CROSS-REFERENCE TO RELATED APPLICATIONS 
5 [0001] This application is related to EP-A-0 993 804, the disclosure of which is incorporated herein by reference. 
FIELD OF THE INVENTION 

[0002] The present invention relates generally to non-contact tracking of objects using a magnetic field, and specif- 
10 ically to counteracting the effect of an intruding field- responsive article in the field. 

BACKGROUND OF THE INVENTION 

[0003] Non-contact electromagnetic tracking systems are well known in the art, with a wide range of applications. 

is [0004] U.S. Patent No. 5,391 ,1 99, whose disclosure is incorporated herein by reference, describes a system for gen- 
erating three-dimensional location information regarding a medical probe or catheter. A sensor coil is placed in the 
catheter and generates signals in response to externally applied magnetic fields. The magnetic fields are generated 
by three radiator coils, fixed to an external reference frame in known, mutually spaced locations. The amplitudes of 
the signals generated in response to each of the radiator coil fields are detected and used to compute the location of 

20 the sensor coil. Each radiator coil is preferably driven by driver circuitry to generate a field at a known frequency, distinct 
from that of other radiator coils, so that the signals generated by the sensor coil may be separated by frequency into 
components corresponding to the different radiator coils. 

[0005] PCT patent publication WO96/05768 whose disclosure is incorporated herein by reference, describes a sys- 
tem that generates six-dimensional position and orientation information regarding the tip of a catheter. This system 

25 uses a plurality of sensor coils adjacent to a locatable site in the catheter, for example near its distal end, and a plurality 
of radiator coils fixed in an external reference frame. These coils generate signals in response to magnetic fields 
generated by the radiator coils, which signals allow for the computation of six location and orientation coordinates. As 
in the case of the '539 patent application described above, the radiator coils preferably operate simultaneously at 
different frequencies, for example at 1000, 2000 and 3000 Hz, respectively. 

30 [0006] The above tracking systems rely on separation of position-responsive signals into components, most typically 
frequency components, wherein each such component is assumed to correspond uniquely to a single radiator coil, in 
a known position, radiating a magnetic field having a regular, well-defined spatial distribution. In practice, however, 
when a metal or other magnetically-responsive article is brought into the vicinity of the catheter or other object being 
tracked, the magnetic fields generated in this vicinity by the radiator coils are distorted. For example, the radiator coil's 

35 magnetic field may generate eddy currents in such an article, and the eddy currents will then cause a parasitic magnetic 
field to be radiated. Such parasitic fields and other types of distortion can lead to errors in determining the position of 
the object being tracked. 

[0007] U. S. Patent 5,767,669 to Hansen et al., whose disclosure is incorporated herein by reference, describes a 
method for subtracting eddy current distortions produced in a magnetic tracking system. The system utilizes pulsed 
40 magnetic fields from a plurality of generators, and the presence of eddy currents is detected by measuring rates of 
change of currents generated in sensor coils used for tracking. The eddy currents are compensated for by adjusting 
the duration of the magnetic pulses. 

[0008] U. S. Patent 4,945,305 to Blood, whose disclosure is incorporated herein by reference, describes a tracking 
system which avoids the problems of eddy currents by using pulsed DC magnetic fields. Sensors which are able to 
45 detect DC fields are used in the system, and eddy currents are detected and adjusted for by utilizing the decay char- 
acteristics and the amplitudes of the eddy currents. 

[0009] EP-A2-0 964 261 , to Dumoulin, whose disclosure is incorporated herein by reference, describes systems for 
compensating for eddy currents in a tracking system using alternating magnetic field generators. In a first system the 
eddy currents are compensated for by first calibrating the system free from eddy currents, and then modifying the fields 
so generated when the eddy currents are detected. In a second system the eddy currents are nullified by using one or 
more shielding coils placed near the generators. 

[0010] Fig. 1 is a graph showing a relation of the permeability u of a ferromagnetic material in a magnetic field vs. 
frequency f at which the field is being generated, as is known in the art. Permeability ^ is a factor in the phase shift 
generated by the magnetic field. The graph applies to a change of the permeability |i of the ferromagnetic material, 
55 generated for an article wherein eddy currents are formed. The change reflects the phase shift in a sensor, caused by 
the article, vs. the frequency f. As is known in the art, additional factors affecting the phase shift are geometry of the 
article, and conductivity of the material. The graph shows a virtually linear change in permeability for small changes 
in frequency. 
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SUMMARY OF THE INVENTION 

[0011] It is an object of some aspects of the present invention to provide methods and apparatus for non-contact 

tracking of an object in an energy field-in the presence of an article which interferes with the field. 

[0012] It is another object of some aspects of the present invention to provide methods and apparatus for minimizing 

the effect of an article which interferes with an energy field used for non-contact tracking of an object 

[0013] In a preferred embodiment of the present invention, an object tracking system comprises one or more sensor 

coils adjacent to a locatable point on an object being tracked, and one or more radiator coils, which generate alternating 

energy fields comprising magnetic fields, in a vicinity of the object when driven by respective alternating electrical 

currents. For each radiator coil, a frequency of its alternating electrical current is scanned through a plurality of values 

so that, at any specific time, each of the radiator coils radiates at a frequency which is different from the frequencies 

with which the other radiator coils are radiating. 

[0014] The sensor coils generate electrical signals responsive to the magnetic fields, which signals are received by 
signal processing circuitry and analyzed by a computer or other processor. When a metal or other field-responsive 
article is in the vicinity of the object, the signals typically include position signal components responsive to the magnetic 
fields generated by the radiator coils at their respective instantaneous driving frequencies, and parasitic signal com- 
ponents responsive to parasitic magnetic fields generated due to the article. The parasitic components are typically 
equal in frequency to the instantaneous frequency of the driving frequency, but are shifted in phase, so that the effect 
at each sensor coil is to produce a combined signal having a phase and an amplitude which are shifted relative to the 
signal when no field-responsive article is present. The phase-shift is a function of the driving frequency, and so will 
vary as each driving frequency is scanned. The computer processes the combined signal to find which frequency 
produces a minimum phase-shift, and thus a minimum effect of the parasitic components, and this frequency is used 
to calculate the position of the object. Varying the driving frequency until the phase shift is a minimum is an effective 
method, not known in the art, for reducing the effect of field-responsive articles on the signal. 

[0015] In preferred embodiments of the present invention, an alternative method is also used in order to find a value 
of the position signal, i.e., of the signal produced without interfering effects of the field-responsive article. Measurements 
of the value of the combined signal are made at a plurality of frequencies. The values obtained are used to solve a 
plurality of simultaneous equations comprising the position signal as one of the unknowns in the equations. Thus, 
varying the driving frequency enables the position signal to be determined in the presence of interfering signals from 
field-responsive articles. 

[0016] The present invention relies on the fact that parasitic magnetic fields, generated by metal or other field-re- 
sponsive articles that receive and re-radiate energy from a radiator coil magnetic field are typically at the same fre- 
quency as the radiator coil field, but are shifted in phase relative thereto. The phase shift and the amplitudes of the 
parasitic fields generally depend on properties of the article, including dielectric constant, magnetic permeability and 
geometrical shape. However, both the phase shift and the amplitude of the parasitic fields can be assumed to be linearly 
dependent on the value of the frequency generating the parasitic field. 

[0017] There is therefore provided, according to a preferred embodiment of the present invention, a method for 
tracking an object including: 

producing an unperturbed energy field at a plurality of predetermined frequencies in the vicinity of the object; 
determining a characteristic of a perturbing energy field induced responsive to the unperturbed field, due to intro- 
duction of an article responsive to the unperturbed field into the vicinity of the object; 

receiving a plurality of resultant signals responsive to the unperturbed and perturbing energy fields generated at 
a location of the object after introduction of the article; 

determining an optimal frequency for the unperturbed energy field from amongst the plurality of predetermined 
frequencies responsive to a parameter of the resultant signals; and 

determining spatial coordinates of the object responsive to the resultant signal at the optimal frequency. 

[001 8] Preferably, producing the unperturbed energy field at the plurality of predetermined frequencies includes scan- 
ning the frequencies sequentially. 

[0019] Further preferably, producing the unperturbed energy field at the plurality of predetermined frequencies in- 
cludes multiplexing at least some of the frequencies. 
[0020] Preferably, receiving the plurality of resultant signals includes: 

measuring a baseline phase value ty^ of each of the plurality of resultant signals at the respective plurality of 
predetermined frequencies before introduction of the article; and 

measuring a phase shift (J) total at the respective plurality of predetermined frequencies after introduction of the 
article, so that the parameter comprises a term \$ Xo al - tyj for each of the plurality of predetermined frequencies; and 
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wherefn determining the optimal frequency includes determining a frequency co at which \§ XoXa * - <|>J is a minimum. 
[0021] Preferably, determining spatial coordinates of the object includes determining spatial coordinates responsive 
to an amplitude of a signal |MJ at the frequency co. 

[0022] Further preferably, determining spatial coordinates of the object includes determining spatial coordinates re- 
sponsive to a phase of a signal at the frequency co. 

[0023] Preferably, producing the energy fields includes producing magnetic fields. Preferably, receiving the signals 
includes receiving electrical signals which are generated responsive to the magnetic fields. 

[0024] There is further provided, according to a preferred embodiment of the present invention, a method for tracking 
an object, including: 



producing an unperturbed energy field comprising a plurality of predetermined frequencies in the vicinity of the 
object; 

producing a perturbing energy field by introduction of an article responsive to the unperturbed field into the vicinity 
of the object; 

15 receiving a respective plurality of signals responsive to the unperturbed and perturbing energy fields generated at 

a location of the object after introduction of the article; and 

determining one or more factors conditional on spatial coordinates of the object responsive to the plurality of signals 
and the respective frequencies. 

20 [0025] Preferably, determining the one or more factors includes: 

assuming a phasor A w of a signal responsive to the unperturbed energy field and a phasor A'^ of a signal respon- 
sive to the perturbing energy field to be directly proportional to a plurality of predetermined currents generating 
the fields; and 

25 assuming a phase $ w of the signal responsive to the unperturbed energy field and a phase tyj of the signal re- 

sponsive to the perturbing energy field to be linearly dependent on the plurality of predetermined frequencies. 

[0026] Preferably, the plurality of frequencies includes at least four frequencies, and the one or more factors include 
the spatial coordinates of the object. 
so [0027] Preferably, receiving the plurality of signals comprises receiving at least four values of a signal Mj at the at 
least four frequencies, and determining the one or more factors includes: 

determining a value of a position signal amplitude Aq, generated responsive to the unperturbed energy field, by 
substituting respective values of the signal M, into an equation 

35 

— - i$>' 

Mj = Aj + afe 

wherein M , is a phasor representing a measured field, A, is a phasor representing the unperturbed field, aj' 
*o represents an amplitude of the perturbing field, represents a phase of the perturbing field, and i represents at 

least four numbers respectively corresponding to the at least four frequencies, so as to generate at least four 
equations; and 

solving the at least four equations for the position signal amplitude Aq. 

45 [0028] There is further provided, according to a preferred embodiment of the present invention, object tracking ap- 
paratus, comprising: 

a radiator, which generates an energy field at a plurality of predetermined frequencies in the vicinity of the object; 
a sensor, fixed to the object, which generates a plurality of signals responsive to the energy field and to an interfering 
50 article responsive to the energy field; and 

signal processing circuitry, which receives the plurality of signals from the sensor and determines an optimal fre- 
quency for the energy field from amongst the plurality of predetermined frequencies responsive to a parameter of 
the signals, and which determines position coordinates of the object responsive to the signal at the optimal fre- 
quency. 



[0029] Preferably, the radiator generates the energy field at the plurality of predetermined frequencies by scanning 
the frequencies sequentially. 

[0030] Further preferably, the radiator generates the energy field at the plurality of predetermined frequencies by 
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multiplexing at least some of the frequencies. 

[0031] Preferably, the parameter includes a phase shift, and the optimal frequency includes the frequency where the 
phase shift is a minimum. 

[0032] Preferably, the signal processing circuitry determines the position coordinates of the object responsive to an 
s amplitude of one of the plurality of signals at the frequency where the phase shift is a minimum. 
[0033] Preferably, the energy field includes a magnetic field. 

[0034] Preferably, the plurality of signals include a plurality of electrical signals which are generated responsive to 
the magnetic field. 

[0035] There is further provided, according to a preferred embodiment of the present invention, object tracking ap- 
10 paratus, including: 

a radiator, which generates an energy field including a plurality of predetermined frequencies in the vicinity of the 
object; 

a sensor, fixed to the object, which generates a respective plurality of signals responsive to the energy field and 
is to an interfering article responsive to the energy field; and 

signal processing circuitry, which receives the plurality of signals from the sensor and determines one or more 
factors conditional on spatial coordinates of the object responsive to the signals and their corresponding frequen- 
cies. 

20 [0036] Preferably, the plurality of frequencies includes at least four frequencies, and wherein the one or more factors 
comprise the spatial coordinates of the object. 

[0037] The present invention will be more fully understood from the following detailed description of the preferred 
embodiments thereof, taken together with the drawings, in which: 

25 BRIEF DESCRIPTION OF THE DRAWINGS 

[0038] 

Fig. 1 is a graph showing a relation of the permeability |x of a ferromagnetic material in a magnetic field vs. frequency 
so f at which the field is being generated, as is known in the art; 

Fig. 2 schematically illustrates a system for tracking a probe, such as a catheter for medical use, according to a 
preferred embodiment of the present invention; 

Fig. 3 is a vector diagram illustrating a relation between position and parasitic components of a signal generated 
in the system of Fig. 2, according to a preferred embodiment of the present invention; and 
35 Fig. 4 is a schematic flow chart showing a method for choosing frequencies to track the probe of the system of 

Fig. 2, according to a preferred embodiment of the present invention. 

DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 

40 [0039] Reference is now made to Fig. 2, which schematically illustrates a system 10 for tracking a probe 20, such 
as a catheter for medical use, according to a preferred embodiment of the present invention. The operation of systems 
generally similar to system 1 0 are described in detail in the above-mentioned EP-A-0 993 804, U.S. Patent 5,391 ,1 99, 
and PCT patent publication WO/96/05768, which are incorporated herein by reference. System 1 0 comprises a plurality 
of radiator coils 22, 24 and 26. These coils generate respective magnetic fields 7? 2 and # 3 , at respective sets of 

<5 frequencies {©J, {co 2 } and {003}, in the vicinity of probe 20. Each frequency set {o^}, {eo 2 } and {to 3 }, comprises a plurality 
of individual frequencies. Most preferably, each of frequencies f ( in sets {co-,}, {o> 2 } and {co 3 } is simply divisible by a 
common frequency f 0 , i.e., fj = k t f 0 where k 1 is a whole number. In this case, a sampling period AT circuitry receiving 
radiated signals, which signals are explained in more detail hereinbelow, is preferably given by the following equation: 




(1) 



where k 2 represents a whole number. 
[0040] Typical values of frequencies fj in sets {o^}, {co 2 } and {co 3 } comprise frequencies in the range 1 00 Hz - 20 kHz. 
[0041] System 10 further comprises substantially similar variable frequency driver circuitry 30, 32 and 33, coupled 
to each of the radiator coils, which drive coils 22, 24 and 26 at the respective sets of frequencies {g^}, {cd 2 } and {0)3}. 
Most preferably, the sets of frequencies {coj, {co 2 } and {©3} at which the coils radiate are set by computer 36. Further 
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most preferably, at least some of the frequencies of each driver circuitry are multiplexed together, and after acquisition 
the resulting signals are analyzed in signal processing circuitry 34, as described in more detail below. Alternatively, 
the frequency of each driver circuitry is varied according to some other method known in the art, such as scanning the 
frequencies cyclically over time, and/or using one or more other methods of time multiplexing. Whichever method is 
used to vary the frequencies, at any instant in time a frequency radiated by a specific coil is set to be different from the 
frequency or frequencies radiated by all the other coils. 

[0042] The probe includes sensor coils 27, 28 and 29, which generate electrical current signals in response to the 
magnetic fields. At any instant in time these signals comprise components of the specific frequencies {c^}, {a> 2 } and 
{0)3} being generated, whose respective amplitudes are dependent on the position and orientation of probe 20. The 
signals generated by sensor coils 27, 28 and 29 are preferably received and processed by signal processing circuitry 
34 and then used by computer 36 to calculate position and orientation coordinates of probe 20. 
[0043] Fig. 2 shows three radiator coils 22, 24 and 26 and three sensor coils 27, 28 and 29 in a probe 20. It will be 
understood, however, that the present invention is equally applicable to tracking systems comprising one, two, four or 
more radiator coils and one, two or more sensor coils. For example, the present invention applies to a single axis 
system comprising one sensor coil, in which case the system most preferably comprises nine radiator coils. 
[0044] In the absence of parasitic effects, the signals generated by sensor coils 27, 28 and 29 at any of frequencies 
{0^} are proportional to the amplitude of the time derivative of the projection of field /?, at probe 20 along the respective 
axes of the sensor coils. The signals generated at any of frequencies {0^} and {0)3} are similarly proportional to the 
projections of Vi 2 and Vi 3 . Parasitic effects that may arise due to mutual inductance among the radiator coils are pref- 
erably substantially eliminated, as disclosed, for example, in PCT patent application no. PCT/IL/00100, whose disclo- 
sure is incorporated herein by reference. 

[0045] Since the direction and amplitude of the magnetic field due to any one of radiator coils 22, 24 and 26 can be 
calculated easily using methods known in the art, the sensor coil signals due to the respective radiator coil field may 
be directly related to the sensor coil's distance from and orientation relative to the radiator coil. It will also be appreciated 
that in the absence of parasitic magnetic fields, such as will be described below, the phase of the signal at each specific 
frequency comprised in {0^ }, {o) 2 } and {0)3} is substantially constant relative to the phase of the magnetic field generated 
by radiator coils 22, 24, 26, and depends on the position and orientation of sensor coils 27, 28, 29. 
[0046] As shown in Fig. 2, however, when a metal or magnetic field-responsive article, for example a surgical tool 
40, is introduced into the vicinity of probe 20, the article will generally receive energy from unperturbed fields ~fi 2 
and H 3f and re-radiate perturbing parasitic magnetic fields, 7?' r , #' 2 and 7?' 3 , at the specific frequencies from sets {o^}, 
{cd 2 } and {0)3} which are being generated. Generally the phases of the parasitic fields will be shifted relative to the 
radiator coil fields by phase angles <|> 1 , ) <|> 2 ' and <f> 3 \ respectively. The phases and amplitudes of the parasitic fields 
generally depend on properties of tool 40, including its dielectric constant, magnetic permeability, geometrical shape 
and orientation relative to the radiator coils. The phases and amplitudes of the parasitic fields are also a function of 
the specific frequencies being generated. 

[0047] Fig. 3 is a vector diagram illustrating a relation between the position and parasitic signal components, for 
radiation from radiator coil 22 at sensor coil 27, according to a preferred embodiment of the present invention. Coil 27 
generates a set of frequencies, responsive to the frequencies generated by radiator coils 22, 24 and 26, which are 
transferred to signal processing circuitry 34. Circuitry 34 separates the received signal into constituent frequencies, 
and recovers the amplitude and phase of each frequency, which are used as described hereinbelow with reference to 
Fig. 3. In the interests of simplicity, unless indicated otherwise the following explanation refers to sensor coil 27, although 
it will be appreciated that sensor coils 28, and 29 behave substantially as coil 27. 

[0048] Signal vector 50, having an amplitude |M {(0} | and a phase fy^ 0131 , represents a signal M {a) j received from 
sensor coil 27 at a set of frequencies {g>}. Vector 50 is the vector sum of position signal component vector 52 and 
parasitic signal component vector 54. Vectors 50, 52, and 54 are referenced in phase to a current l {ft) } in coil 27. Position 
signal component 52 has amplitude A {(o} and a substantially constant baseline phase at frequency {co}. At frequency 
{go} parasitic signal component 54 has a phase shift from the baseline of $ {ai} ' and an amplitude A {(0) \ Unless indicated 
otherwise, the following explanation considers one specific frequency, herein termed co, although it will be appreciated 
that the explanation applies to all frequencies generated in coil 27. 

[0049] The total combined signal M w received from sensor coil 27, including both position and parasitic signal com- 
ponents may generally be expressed as: 



Moo = Aco + A'^ (2a) 

so that the amplitude |MJ is given by 



6 



EP 1 203 560 A2 



10 



15 



IMJ-IA^ + A'J (2b) 

wherein A w and A*^ are the phasors of the position signal component and the parasitic signal component re- 
spectively, at frequency co. 

[0050] It will be observed in equation (2b) and from Fig. 3 that for each of the signal frequency components |iVy, 
the superposition of the parasitic signal component will cause a phase shift in the total detected signal, relative to the 
signal phase in the absence of metal tool 40, given by: 



<})£ tal = arctan 



A(psin 4>q) + Aq) sin 
Ato cos <))(„) + A W cos ^ 



(3) 



[0051] In preferred embodiments of the present invention, signal processing circuitry 34 and computer 36 detect and 
record baseline phases <t> w for all different frequencies received from sensor coils 27, 28 and 29, or for other systems 
described hereinabove such as the single axis system, in the absence of any metal or other interfering magnetic field- 
20 responsive objects in the vicinity of probe 20. Alternatively, undisturbed phases of the position signal components may 
have been determined in advance for system 1 0 or are known based on the operation of the system. When metal tool 
40 is introduced into the vicinity of probe 20, the phase shift due to the parasitic components engendered thereby in 
the signals is measured at each separate frequency. 

[0052] Fig. 4 is a schematic flow chart showing a method for choosing frequencies to track probe 20, according to 

25 a preferred embodiment of the present invention. For simplicity, the following description refers only to radiator coil 22 
and sensor coil 27, but it will be understood that the method shown in Fig. 4 applies to any combination of a radiator 
coil and a sensor coil in system 10. In an initialization phase a baseline phase value ty^ is measured by circuitry 34 
from the voltages and currents induced in coil 27. The baseline phase value ^ is measured at each of the frequencies 
of {a^} and each value is recorded in computer 36. During an operation phase the value of <|> total and the absolute value 

30 of the difference, |<|> total - <|>J t is measured and recorded for each frequency of In trie event that the absolute 
difference is not equal to zero, indicating that a parasitic signal component due to tool 40 is present, computer 36 
selects the frequency having the smallest absolute difference. This frequency is used when evaluating |M J in equation 
(2b). It will be appreciated that applying the method described hereinabove to all combinations of radiator coils and 
sensor coils in system 10 enables a complete determination of the position and orientation of probe 20. 

35 [0053] As stated hereinabove, equations (2a) and (2b) apply for frequency co applied to sensor coil 27. As the fre- 
quency co is varied, values of A^, <t> w , A^' and tyj varv - As is known in the art, values of A w and A^' are directly proportional 
to the current at which the specific radiator coil, assumed herein to be radiator coil 22, generating the field is being 
driven and which sensor coil 27 is detecting. Thus = (SAq, and A^' = pA 0 ' where p is a constant, oo 0 is an arbitrary 
frequency in {coj}, and Aq and A© are the amplitudes of the position and parasitic signal components at frequency coq. 

40 Also, for small variations of frequency co, tyj is linearly dependent on frequency co, so that 



45 where Aco = co - co 0 , 



y is a constant, corresponding to a value of the derivative 



6^ is the parasitic phase and 



A© 



$0 is the position phase. 
50 At a particular frequency coj, equation (2a) can be rewritten as: 

Mj = Aj + A/ (5a) 

55 where Mj is a phasor representing the measured field at C0j, 

Aj is a phasor representing the unperturbed field, and 

Aj' is a phasor representing the perturbing field due to tool 40. Equation (5a) can be rewritten as follows: 
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Mi = Aj + ai'e'*' (5b) 

where a- and are'the "perturbing amplitude arid phase at 6^. " " " 

Equation (5b) can also be rewritten: 

M ( = Pi |A 0 |e l V Pi a 0 -e^ + ^-'° o) . (5c) 

where a^ and $q are the perturbing amplitude and phase, and Aq is the unperturbed amplitude, at cqq, 
<(>i is the unperturbed phase shift at coi, 



ft* 



M 0 



and 7 = — . 

20 9(0 



[0054] In equation (5c) A 0 , ao*, and <j> 0 ' are unknown, and ty v p j( y, o^, and ©j are known, or in the case of y may be 
found from one other frequency apart from co 0 and o)j by_using the graph of Fig. 1 relating permeability to frequency. 
Alternatively, y may be assumed to be unknown. Thus, if M j is measured at four known separate frequencies, equation 
(5c) can be solved for Aq, the position signal component Most preferably, frequencies {a^}, {c^} and {©3} in system 
10 comprise more than four separate frequencies, so that a plurality of values of Aq can be determined, and a final 
value of Aq calculated by one of the processes of averaging known in the art. Alternatively or additionally, when fre- 
quencies {g^}, {o) 2 } and {g^} in system 1 0 comprise more than four separate frequencies, equation (5c) may be adapted 
to comprise other parameters describing at least some Aq, ao', and ty 0 '. For example, values of Aq and ao' can be 
assumed to depend on frequency co in a linear or a nonlinear manner, and appropriate constants can be included in 
equation (5c), as is known in the art. 

[0055] It will thus be appreciated that by varying the excitation frequency of each radiator coil, and measuring the 
total signal generated in each sensor coil at these frequencies, the position component of the signal can be determined 
regardless of the presence of parasitic components. It will also be appreciated that varying the excitation frequency of 
each radiator coil by a plurality of frequencies, wherein the plurality is fewer than four, will give useful information 
regarding factors associated with tracking objects in the presence of interfering articles. 

[0056] It will further be appreciated that the preferred embodiments described above are cited by way of example, 
and that the present invention is not limited to what has been particularly shown and described hereinabove. Rather, 
the scope of the present invention includes both combinations and subcombinations of the various features described 
hereinabove, as well as variations and modifications thereof which would occur to persons skilled in the art upon reading 
the foregoing description and which are not disclosed in the prior art. 



Claims 

1. A method for tracking an object, comprising: 



producing an unperturbed energy field at a plurality of predetermined frequencies in the vicinity of the object; 
determining a characteristic of a perturbing energy field induced responsive to the unperturbed field, due to 
introduction of an article responsive to the unperturbed field into the vicinity of the object; 
receiving a plurality of resultant signals responsive to the unperturbed and perturbing energy fields generated 
at a location of the object after introduction of the article; 

determining an optimal frequency for the unperturbed energy field from amongst the plurality of predetermined 
frequencies responsive to a parameter of the resultant signals; and 

determining spatial coordinates of the object responsive to the resultant signal at the optimal frequency. 

2. A method according to claim 1 , wherein producing the unperturbed energy field at the plurality of predetermined 
frequencies comprises scanning the frequencies sequentially. 
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3. A method according to claim 1 , wherein producing the unperturbed energy field at the plurality of predetermined 
frequencies comprises multiplexing at least some of the frequencies. 

4. A method according to claim 1 , wherein receiving the plurality of resultant signals comprises: 

measuring a baseline phase value ^ of each of the plurality of resultant signals at the respective plurality of 
predetermined frequencies before introduction of the article; and 

measuring a phase shift (|> total at the respective plurality of predetermined frequencies after introduction of the 
article, so that the parameter comprises a term |(t> total - $ J for each of the plurality of predetermined frequencies; 
and w 

wherein determining the optimal frequency comprises determining a frequency o> at which |<J) total - <|>J is a 
minimum. tL 

5. A method according to claim 4, wherein determining spatial coordinates of the object comprises determining spatial 
coordinates responsive to an amplitude of a signal |MJ at the frequency co. 

6. A method according to claim 4, wherein determining spatial coordinates of the object comprises determining spatial 
coordinates responsive to a phase of a signal at the frequency co. 

7. A method according to claim 1 , wherein producing the energy fields comprises producing magnetic fields. 

8. A method according to claim 7, wherein receiving the signals comprises receiving electrical signals which are 
generated responsive to the magnetic fields. 

9. A method for tracking an object, comprising: 

producing an unperturbed energy field comprising a plurality of predetermined frequencies in the vicinity of 
the object; 

producing a perturbing energy field by introduction of an article responsive to the unperturbed field into the 
vicinity of the object; 

receiving a respective plurality of signals responsive to the unperturbed and perturbing energy fields generated 
at a location of the object after introduction of the article; and 

determining one or more factors conditional on spatial coordinates of the object responsive to the plurality of 
signals and the respective frequencies. 

10. A method according to claim 9, wherein determining the one or more factors comprises: 

assuming a phasor A w of a signal responsive to the unperturbed energy field and a phasor A'^ of a signal 
responsive to the perturbing energy field to be directly proportional to a plurality of predetermined currents 
generating the fields; and 

assuming a phase <t> w of the signal responsive to the unperturbed energy field and a phase ^ of the signal 
responsive to the perturbing energy f ield to be linearly dependent on the plurality of predetermined f requencies. 

1 1 . A method according to claim 9, wherein the plurality of frequencies comprises at least four frequencies, and wherein 
the one or more factors comprise the spatial coordinates of the object. 

12. A method according to claim 11 , wherein receiving the plurality of signals comprises receiving at least four values 
of a signal Mj at the at least four frequencies, and wherein determining the one or more factors comprises: 

determining a value of a position signal amplitude Aq, generated responsive to the unperturbed energy field, 
by substituting respective values of the signal Mj into an equation 

— - ify* 
Mj = Aj + aj'e 

wherein Mj is a phasor representing a measured field, Aj is a phasor representing the unperturbed field, a{ 
represents an amplitude of the perturbing field, <fc' represents a phase of the perturbing field, and i represents 
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at least four numbers respectively corresponding to the at least four frequencies, so as to generate at least 
four equations; and 

solving the at least four equations for the position signal amplitude Aq. 

13. Object tracking apparatus, comprising: 

a radiator, which generates an energy field at a plurality of predetermined frequencies in the vicinity of the 
object; 

a sensor, fixed to the object, which generates a plurality of signals responsive to the energy field and to an 
interfering article responsive to the energy field; and 

signal processing circuitry, which receives the plurality of signals from the sensor and determines an optimal 
frequency for the energy field from amongst the plurality of predetermined frequencies responsive to a pa- 
rameter of the signals, and which determines position coordinates of the object responsive to the signal at the 
optimal frequency. 

14. Apparatus according to claim 13, wherein the radiator generates the energy field at the plurality of predetermined 
frequencies by scanning the frequencies sequentially. 

15. Apparatus according to claim 13, wherein the radiator generates the energy field at the plurality of predetermined 
frequencies by multiplexing at least some of the frequencies. 

16. Apparatus according to claim 13, wherein the parameter comprises a phase shift, and wherein the optimal fre- 
quency comprises the frequency where the phase shift is a minimum. 

17. Apparatus according to claim 16, wherein the signal processing circuitry determines the position coordinates of 
the object responsive to an amplitude of one of the plurality of signals at the frequency where the phase shift is a 
minimum. 

18. Apparatus according to claim 13, wherein the energy field comprises a magnetic field. 

19. Apparatus according to claim 18, wherein the plurality of signals comprise a plurality of electrical signals which 
are generated responsive to the magnetic field. 

20. Object tracking apparatus, comprising: 

a radiator, which generates an energy field comprising a plurality of predetermined frequencies in the vicinity 
of the object; 

a sensor, fixed to the object, which generates a respective plurality of signals responsive to the energy field 
and to an interfering article responsive to the energy field; and 

signal processing circuitry, which receives the plurality of signals from the sensor and determines one or more 
factors conditional on spatial coordinates of the object responsive to the signals and their corresponding fre- 
quencies. 

21. Apparatus according to claim 20, wherein the plurality of frequencies comprises at feast four frequencies, and 
wherein the one or more factors comprise the spatial coordinates of the object. 
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